Current collapse in tunneling transport through benzene 
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We investigate the electrical transport through a system of benzene coupled to metal electrodes 
by electron tunneling. Using electronic structure calculations, a semi-quantitative model for the 
TT electrons of the benzene is derived that includes general two-body interactions. After exact 
diagonalization of the benzene model the transport is computed using perturbation theory for weak 
electrode-benzene coupling (golden rule approximation). We include the effect of an applied electric 
field on the molecular states, as well as radiative relaxation. We predict a current collapse and 
strong negative differential conductance due to a "blocking" state when the electrode is coupled to 
the para-position of benzene. In contrast, for coupling to the meta-position, a series of steps in the 
I ~ V curve is found. 

PACS numbers: 73.63.-b, 73.23.Hk, 73.22.-f 



Introduction. Single molecule electronic devices offer 
exciting perspectives for further minituarization of elec- 
tronic circuits with a potentially large impact in applica- 
tions. Several experiments have demonstrated the pos- 
sibility to attach individual molecules to leads and to 
measure the electrical transport. |Q, ||, ^, ^ In contrast 
to single electron transistors (SETs) based on quantum 
dots 1^ the electronic structure of molecular devices can 
be chemically 'designed' for specific applications. When 
the molecule is coupled weakly to the electrodes, i.e. via 
electron tunneling, charging effects, semi-classically de- 
termined by the small capacitance of the molecule, be- 
come important. The interplay of charging effects with 
the specific structure of the molecular orbitals leads to 
nontrivial current voltage (/ — V) characteristics Q . 
Very recent experiments demonstrated both Coulomb 
blockade and the Kondo effect in three terminal trans- 
port through a single molecular level j|, jo) . 

Using benzene as a prototypical example we investi- 
gate novel effects that arise when transport through sev- 
eral competing electronic configurations becomes possi- 
ble. We derive a semi-quantitative model for the con- 
ducting many-body states of the system from electronic 
structure calculations. For weak coupling to the electrode 
we compute transport within the golden rule approxima- 
tion (sequential tunneling) and include screening of the 
applied electric field as well as radiative transitions be- 
tween the electronic states of the molecule. We predict a 
current collapse in the I — V curve and strong negative 
differential conductance (NDC) due to the occurrence of 
a "blocking" state when the molecule is coupled to the 
electrode at the para-position. For coupling at the meta- 
position, the I — V curve displays a series of steps, but 
no NDC. We demonstrate how the specific spatial struc- 
ture of the molecular orbitals qualitatively determines 
electronic transport. Finally, we discuss the limits of the 
model and the impact of disorder and symmetry breaking 
effects likely encountered in an experimental realization. 



The Model. To perform transport calculations in the 
weakly coupled regime, we extract an effective model 
from electronic structure calculations of the molecule. 
For benzene transport is assumed to be dominated by 
the TT-electron system, but generalizations are straight- 
forward. We first perform Hartree-Fock calculations in 
a suitable basis and then transform the Hamiltonian to 
the molecular orbital basis. We then integrate out the a 
electrons of the system arriving at an effective interacting 
model Hamiltonian for the tt electrons of the system in 
the presence of the atomic cores and the "frozen" density 
of the a electrons: 
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The second quantized operators clj^,Ca^i create/destroy 
electrons of spin a in orthogonalized Wannier-like or- 
bitals <i>i centered at the carbon atoms. While this 
model neglects ct-tt mixing for certain excited states of 
the molecule 11 , its parameters account for the de- 
tailed electronic structure of the molecule. For the cur- 
rent work, we compute the model parameters only once, 
using an augmented double-zeta quality atomic natu- 
ral orbital (ANO) basis set that was truncated to con- 
tain only one 2p-shell for each carbon. When applying 
a bias over the molecule this neglects the higher order 
effect of field screening by a electrons and its impact 
on the TT electrons. Note that in our approach there is 
no over-counting of interaction terms like in "interaction 
enhanced" density functional treatments in solid state 
physics ll|,|l|. 

We find that the low-energy spectrum of benzene ob- 
tained from the diagonalization of the effective Hamil- 
ton operator eq. |l| compares favorably with the spectrum 
directly obtained from accurate multi-reference configu- 
ration interaction calculations ||ll|, |l^. The remaining 
differences can be understood by the lack of a-ir mixing 
and do not qualitatively affect the transport properties. 
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The restriction to a pure 7r-electron system is more severe 
for the charged states, as discussed in some detail below. 

To account for the effect of an external bias potential 
yext ^Yie electrons we include a term 

Ht^as = e ^ V,f'cl,c,,, (2) 

with V^f * = / dr$i(r)F'==^*(f)$j (r) and F^^*(f) {Vl + 
VR)/2 — Vbias{^/L) in the Hamiltonian. Vl,r is the chem- 
ical potential in the left/right electrode, distances are 
measured from the center of the molecule. The ma- 
trix elements K';^* also result from the electronic struc- 
ture calculation described above. We consider a bias 
Vbias = Vl — Vr aligned with the transport direction 
(x-direction) . The parameter L indicates the length over 
which the external bias falls off (we choose L = OAnm). 

Diagonalization and Transport Calculation. To com- 
pute the transport properties of the system in the weak 
coupling limit we diagonalize the Hamiltonian eq. |^ in 
the appropriate charge/spin/symmetry sectors. For the 
current model, we have one "single-particle" state (Wan- 
nier state ^i) per carbon site i, resulting in a total of 
4^ = 4096 (with iV = 6 sites) many-body states \s). 
For the present model, charge, spin and symmetry adap- 
tation reduce the size of effective Hamiltonian matrices 
sufficiently to permit their diagonalization with standard 
linear algebra packages. Care is needed to handle degen- 
eracies properly. The energy degenerate states need to 
be additionally diagonalized to obtain states with inte- 
ger (half-integer) total spin for states with even (odd) 
number of electrons. This also helps for the transport 
calculation, since the transition matrix elements can be 
summed over the magnetic quantum numbers, i.e. we 
only need to consider one representative of the degener- 
ate spaces (we do not consider application of a magnetic 
field here). 

After diagonalization of the Hamiltonian we have the 
many-body eigenstates \s) with the corresponding ener- 
gies Eg and their total spin Ss ■ We use a Master equation 
approach Q for the occupation probabilities Pg in a sta- 
tionary state. The transition rates T,ss' from state s' to 
s are computed in perturbation theory using the golden 
rule. The "perturbation" is the coupling of the molecule 
to the leads 

Hmol-leads = (^^)^^^ ^ {c]„a]^aa + h.c)j , (3) 

and (optionally) the coupling to electromagnetic fields 
(photons), r is the coupling strength (in units of en- 
ergy) of leads to the benzene and pe is the density of 
states of the electrons in the electrode (assumed con- 
stant). The operators flko-a and their hermitian conju- 
gates destroy/create electrons with momentum k and 
spin a in electrode a = left /right. For simplicity, we 



assume that tunneling is only possible through two "con- 
tact" carbon atoms which we choose to be at the 1 and 
4 (para) positions unless noted otherwise. 

As we do not consider the leads microscopically, the 
coupling of molecule states \s) is determined by the over- 
all coupling strength F and the relative wave function am- 
plitude of the state \s) at the coupling carbon site /. For 
the transition rates we have Sss' = (X]q p=± ^^f' ) + s' 
where E"!! is the tunneling rate to/from electrode a for 
creation (p = -|-) or destruction {p = — ) of an electron 
on the molecule. We have 

= TUEs ~ E,) E I E < ^l^ll^' > I' ' (4) 

and a corresponding equation for by replacing fa — > 
1 — /a, where fa is the Fermi function. When includ- 
ing relaxation by radiative transitions, we use the dipole 
approximation with dipole transition moments d^j = 
/ (ir$j(r)r<i>j(r) obtained from the electronic structure 
calculations. The corresponding transition rates are 

^is' = ^i^s - Es.fNb{Es -Es,)\< s\d\s' > p 

where 7Vf,(i?) denotes the equilibrium Bose function. 
Note that for emission Eg — Egt is negative, and 

mHE\) = -ii + mi\E\))- 

The total transition matrix S^s' consists of blocks con- 
necting N and A*" ± 1 electron states (tunneling pro- 
cesses) and blocks from the radiative transitions that do 
not change the electron number on the molecule. Tak- 
ing only one member of the subspace of spin and en- 
ergy degenerate states into account, the rank r of the 
transition matrix is 1716. The stationarity condition 
Ps — can be written as Ags'Ps' — with the matrix 
Ass' = Sss' — X^s" ^s"sSss'- This implies that Ass' has 
an eigenvector with zero eigenvalue, which is the wanted 
solution for Pg. Rather than computing this eigenvector 
by brute force, to speed up the calculation we make use 
of X^ = 1 to eliminate one row/column, thus refor- 
mulating the eigenvector problem into one of solving an 
inhomogeneous linear system of rank r — 1. Eventually, 
the current in the left and right electrode is calculated 
via/„ = eE...'(S:;n'-E:,;P,). ^ 

In this paper, we consider symmetric bias only. Apply- 
ing an asymmetric bias would effectively imply different 
coupling to the left and right leads and has been previ- 
ously discussed [Q. Finally, we note that we must have 
F <C kbT w AOmeV for perturbation theory to apply at 
room temperature T. 

Results: To elucidate the impact of various effects on 
the current we have performed transport calculations 
with and without radiative transitions (relaxation) and 
with and without effect of the applied bias. 

On the left panel of Fig. |l| we present the I — V 
curve obtained without the effect of the applied field 
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FIG. 1: I ~ V characteristics for symmetric bias. Left Panel: 
I —V curve without (solid line) and with radiative relaxation 
(dashed line) without inclusion of bias effect. The current col- 
lapses at a bias when the antisymmetric anion state can be- 
come occupied by radiative relaxation from an excited anion 
state. For coupling of the right electrode at the meta position 
(dash-dotted line) the "blocking state" has finite tunnel cou- 
pling to both electrodes and no current collapse is observed. 
Right Panel: With inclusion of the bias effect (solid line) the 
onset of current is generally shifted to lower bias, but the 
current collapse remains and additional weak NDC occurs at 
larger bias. 



(zero bias electronic structure). Without radiative re- 
laxation (solid line) the I — V curve consists of a series of 
steps of which only a few are well resolved on this scale. 
The first step is associated with the population of the first 
TT* orbital of molecule (molecular charge = — e), an elec- 
tron hops onto the lowest available level and then hops 
off again. At slightly larger bias, a transition of the anion 
to the first excited state of the neutral molecule becomes 
possible, resulting in a slight increase of the current. If 
the bias is sufficiently large this excited state may now 
accept another electron to populate higher excited states 
of the anion or low-lying states of the di-anion, result- 
ing in a rapidly growing cascade of transitions between 
literally hundreds of states of the system. In the model 
considered here the growth of this cascade leads to quasi- 
ohmic behavior above 3.6 eV. In our calculation the first 
states of the di-anion become occupied at about 4.5 eV. 

The inclusion of radiative transitions has a dramatic 
effect on the I — V curve (dashed line). We observe a 
collapse of the current over a substantial range of the 
applied bias (2.1-3.4 eV). The reason for this collapse is 
the population of a "blocking" state in the cascade of 
transitions that becomes possible when exited states of 
the neutral molecule and anion become accessible. Above 
approximately 2.1 eV bias an excited state of the anion 
at about 5.6 eV (see Fig. ||), becomes partially populated 
in the transport cascade. This state can decay by pho- 
ton emission to either a symmetric or an antisymmetric 
many body state (with respect to the plane through the 



3.4 eV 




Neutral benzene Charged benzene anion 

FIG. 2: Sketch of the energetics and symmetry of the rele- 
vant neutral and anion states. At a left bias of 2.1 eV the 
antisymmetric blocking state becomes occupied via radiative 
relaxation. At 3.4 eV electrons can escape the blocking state 
via tunneling to the antisymmetric state of neutral benzene. 



transport axis and perpendicular to the molecular plane) 
of the anion. In the bias range of the current collapse 
this state cannot decay by coupling to the leads, because 
the lowest neutral states are symmetric (see Fig. ||) and 
the tunneling preserves the symmetry. Since there are no 
further radiative transitions possible on the molecule, the 
rate equations contain no draining term from this state. 
As a result, in the stationary state, the probability of 
occupying the "blocking" state is unity and the current 
ceases to fiow. At a larger bias (3.4 eV), the first escape 
channel opens, and the system can decay to the first an- 
tisymmetric excited state of the neutral molecule, which 
can then decay further by photon emission. 

The above calculations simplify the illustration of the 
mechanism of NDC, because the energies of the partici- 
pating states are independent of the applied bias. In the 
physical system the various states couple differently to 
the applied field. The solid line in the right panel of Fig. |l| 
shows the I — V curve of the benzene with radiative re- 
laxation and with the effect of the external potential ap- 
plied parallel to the transport axis. We note a shift of the 
first step in the I — V curve, which is due to differential 
screening effects between the ground state of the neutral 
molecule and the anion. At larger bias, the differential 
effects on all of the states in the cascade result in a sig- 
nificant renormalization of the I — V curve. However, the 
current collapse still occurs, although the voltage window 
is reduced. The onset of quasi-ohmic behaviors occurs at 
higher voltage and a series of weak NDC effects occurs 
due to redistribution of occupation probabilities in favor 
of states with slightly smaller "transmission" . 

If we couple the right electrode to the meta position of 
benzene (left panel, dash-dotted curve) we find no cur- 
rent collapse with or without radiative relaxation or bias 
effect. This can be readily understood by the fact that 
at the meta position the wave function of the formerly 
blocking state is non-vanishing and electrons can tunnel 
out to the right electrode. Consequently, we observe a 
series of current steps similar (but not identical) to the 
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case of coupling to the para position without relaxation. 

Discussion: In the above we have described trans- 
port in an idealized, though semi-quantitative, 7r-electron 
model of weakly coupled benzene and found a dramatic 
suppression of the current in a finite voltage window. The 
mechanism for this effect is the occupation of a so-called 
blocking state of the molecule, which cannot decay for 
symmetry reasons. In the following we will summarily 
discuss its stability under the influence of a number of 
effects not considered in our calculation, but likely to 
be present in a experimental realization. Details will be 
discussed in a forthcoming publication |]l5| . 

(i) Transport through the a-system: The ionization 
potential of benzene is rather high, for this reason no 
cationic states (with either cr or tt holes) are relevant 
for the bias range considered here and attainable in ex- 
periment. However, for the anion, low lying Rydberg 
states are predicted in quantum chemical calculations 
with augmented basis sets containing diffuse functions. 
In these states the additional particle occupies an elec- 
tron cloud of cr symmetry that is smeared out over the 
entire molecule and has significant charge density in the 
center of the ring. Including these states, we expect the 
onset of current in the I — V curve, i.e. the first current 
step, will be due to tunneling into the lowest Rydberg 
state, at about 1 eV. Nevertheless we believe that the 
NDC effect will be preserved. The current collapse is due 
to the competition between a slow process (cascade fill- 
ing the blocking level) with an even slower one (decay of 
the blocking level). NDC occurs when probability den- 
sity is redirected from "conducting" many-body states 
with large tunnel coupling to non-conducting many-body 
states even as the voltage increases. Note that the exis- 
tence of a single "blocking path" is sufficient for NDC, ir- 
respective of the number of competing non-blocking cas- 
cades. 

The addition of Rydberg states below the blocking 
TT state of the anion permit its decay to symmetric or 
antisymmetric Rydberg states of the anion. Note that 
the spectrum of the neutral molecule is unaffected — 
no additional low-lying configurations must be consid- 
ered. Therefore a new radiative cascade on the anion ul- 
timately leads to a symmetric Rydberg state, which can 
couple to the ground state configuration of the neutral 
molecule. The amplitude for this relaxation, however, 
is slow, because the transition dipole moment from the 
blocking (tt) configuration to the Rydberg (cr) configura- 
tion is small jist (involving a dipole transition moment 
perpendicular to the benzene plane). 

(ii) Symmetry breaking: In most experimental realiza- 
tions for two terminal transport, with the possible excep- 
tion of STM contacts ||l^ , the alignment of the molecule 
with respect to the electrodes is uncontrolled. It is there- 
fore likely that the applied bias is not aligned with the 
transport direction. In this case the degenerate states 
of the anion mix, and transport is no longer completely 



blocked. However, according to the rate argument pre- 
sented above, it is still impeded, and a reduced NDC 
remains observable. The same consideration holds true 
for slow vibrational modes of the molecule, like breathing 
and buckling jl^ . Fast vibrations can provide radiation- 
less relaxation and should be considered at the quantum 
level. They can be particularly important if they gener- 
ate transition matrix elements that allow for dipole for- 
bidden electronic transitions in the undistorted system. 

Summary: In weakly coupled molecules the concept of 
the "blocking" state suggests new opportunities for the 
design of single molecule electronic devices. In this in- 
vestigation we have demonstrated how transport through 
several competing many-body states can result in strong 
NDC and current collapse over a significant voltage win- 
dow. Our calculations predict the existence of this effect 
for a realistic molecule. Note that for Vbias in the vicin- 
ity of the current peak, the application of a gate voltage 
Vg may induce a large asymmetric change in the current, 
which would first rise and then drop again upon scan- 
ning of Vg through a certain range. This behavior could 
be useful in various electronic devices. 
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